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The twin-peak high-frequency quasi-periodic oscillations (HF QPOs), observed in the power spec¬ 
tra of low-mass X-ray binaries, might carry relevant clues about the physics laws reigning close to a 
compact object. Their frequencies are typical of the orbital motion time-scales a few gravitational 
radii away from the compact object. The aim of the manuscript is to propose an intuitive model 
explaining that the energy carried by the lower HF QPO can be related to differences of potential 
energy released by clumps of plasma spiraling in a curved space-time. Our model provides estimates 
on both the size of clumps of matter that can survive to the strong tidal force and energy loaded 
by tides on the clump. We have also obtained some constraints on the mechanical properties of 
the plasma orbiting into the accretion disk. We note that the systematic behavior of the emitted 
energy as function of the central frequency of the lower HF QPO, observed in several sources with a 
neutron star, might give clues related to an innermost stable bound orbit predicted by the General 
Relativity theory in strong field regime. 

PACS numbers: 95.30.Sf, 97.80. Jp, 97.10.Gz, 91.60.Ba 


I. INTRODUCTION 

Since their discovery jlj, the twin-peak high-frequency 
quasi-periodic oscillations (HF QPOs) observed in the 
power spectra of low-mass X-ray binaries (LMXBs) have 
attracted much attention. The HF QPOs have frequen¬ 
cies that correspond to the time-scales of the orbital mo¬ 
tion close to the compact object (milliseconds). They 
are believed to be produced by orbiting inhomogeneities 
of the density of the accretion disk, often referred as to 
clumps/blobs of matter JH3 , or to oscillation modes of 
the accretion disk [^4lo|. Twin-peak HF QPOs might 
be potential probes to disclose the imprints produced 
by the orbiting matter in a strongly curved space-time, 
e.g. the signature of an innermost stable circular orbit 
(ISCO) predicted by General Relativity (GR) in strong 
field [n|, [Hi], as we ll as the modulations of the X-ray 
flux by the precession of the orbits (see Refs. @,1 for 
review). HF QPOs can also be used to constrain fun¬ 
damental quantities of the compact object, such as its 
mass m and angular momentum a 0. Attempts to esti¬ 
mates such quantities in systems with a neutron star (NS 
LMXBs) were done soon after the discovery of HF QPOs 
BUI- Most recently, Motta et al. |Ijl B] gave mea¬ 
surements of both m and a of black hole (BH) LMXBs. 
The observed frequencies of the HF QPOs were linked 
to the relativistic frequencies as in the relativistic pre¬ 
cession model (RPM; Ref. 1), where the upper peak in 
frequency of the twin-peak HF QPOs is linked to the 
Keplerian frequency Vk of the orbiting matter, while the 
lower peak to the periastron precession frequency v p of 
the orbit. Different interpretatio ns |l5l and models link¬ 
ing Vk to the lower HF QPO fl8l - 120j were also proposed. 


* claudio.germana@gmail.com claudio.germana@ufma.br 


The unprecedented possibility of constraining the 
space-time around a compact object through HF QPOs 
and the idea of preferred orbital radii where HF QPOs 
would be produced have stimulated several works based 
on resonance mechanisms. Resonance between relativis¬ 
tic epicyclic and Keplerian frequency of the orbiting mat¬ 
ter has been proposed to justify the observed 3:2 ratio of 
the twin-peak HF QPOs j21H24| . Other resonance mech¬ 
anisms propose forced oscillations in the accretion disk 
induced by the magnetic field of the neutron star [25j . 
or resonance because of coupling between the oscillation 
modes of the orbiting matter in the disk and the spin of 
the neutron star |26l - l30j ]. Also, substantial work was done 
to test the models with the observed data and constrain 
both m and a of the compact object [3ll - f33| . 


Beside being potential probes to study the motion of 
matter in a strongly curved space-time, works were de¬ 
voted to illustrate other properties of the HF QPOs. de 
Avellar et al. [34j , Barret [35] reported time-lags between 
photons from twin-peaks as collected in different energy 
bands in NS LMXBs. Such lags were previously also 
noted in Ref. (3^. Thanks to the larger sample of col¬ 
lected data, de Avellar et al. 34j], Barret HH were able to 
get finer results. These information may shine light on 
the emission mechanisms in such extreme environment 


371. 


Barret et al. [38l. [30l|. Mendez j40[ analyzed the coher¬ 
ence Q of the twin-peaks, defined as Q = i//Azq where v 
is the central frequency of the peak and Azz its width. In 
several NS LMXBs the Q-factor of the lower HF QPO 
displays a steep increase and then an abrupt drop as 
function of the central frequency v of the peak. It was 
proposed that the abrupt drop of the Q-factor might be 
caused by the approach of the oscillation to ISCO [38}, 
thus a geometry-related effect. However, this interpreta¬ 
tion was debated in Refs. [13 l4lj| . If confirmed, a sig¬ 
nature of ISCO is of relevant importance, because it is 
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predicted by GR in strong field regime. In Refs. [38l — 
Eol . E3 | is reported the amplitude of the twin-peak HF 
QPOs as function of their central frequency. Like the 
quality factor Q the behavior of the amplitude (the en¬ 
ergy carried by the QPOs) of the lower HF QPO displays 
a systematic trend with an increase and then a slightly 
steeper decrease for increasing central frequency of the 
peak. The amount of energy released is up to ~ 20 % 
the total flux of the source in NS LMXBs (10 34 — 10 36 
erg/s, from Atoll to Z-sources 0 )- But, how this huge 
quantity of energy is produced is not yet understood. 

Matter orbiting close to a compact object experiences 
a strong tidal force. The effects of the encounter of stars 
with massive black holes were studied by a certain num¬ 
ber of authors [43l - [47j and are of interest for the flaring 
activity observed at the center of galaxies [3 . Reis et al. 
[4^ | reported a QPO in the X-ray flux coming from the 
tidal disruption of a star by a supermassive black hole. 
Most recently, Del Santo et al. [501 linked to the tidal 
disruption of a planet by a white dwarf the X-ray flaring 
activity observed in the source IGR J17361-4441. Subse¬ 
quent analysis revealed a QPO in the X-ray flux produced 
by such event [51] . 

The work done by tides is a significant source of energy al¬ 
ready in our solar system. The volcanism in the Jupiter’s 
moon Io is sustained by tides. Thus, it could be relevant 
to take into account the gravitational energy extracted 
through tides by the orbiting matter close to the com¬ 
pact object in LMXBs. Kostic et al. [52] modeled the 
signal from the tidal disruption of small satellites orbiting 
a Schwarzschild black hole. Germana et al. [ 53 } showed a 
simulated power spectrum with such code. For the first 
time is reported a power spectrum alike to the observed 
ones in LMXBs, both the power law and twin-peaks are 
reproduced. 

The shape of the amplitude of the lower HF QPO, 
as function of its central frequenc y (fi g. 1 in Ref. [42| . 
fig. 3 in Ref. Q and fig. 2 in Ref. [40]), is typical of the 
potential energy released by matter spiraling in a curved 
space-time. We believe that it might be worthwhile to in¬ 
vestigate such issue. In the following we concentrate on 
the lower HF QPO. We show for the first time that the 
energy carried by the lower HF QPO can be accounted for 
by the gravitational energy extracted by spiraling clumps 
of matter undergoing the work done by the strong tidal 
force. We make use of the Schwarzschild potential as 
first-order approach to the problem and check our cal¬ 
culations also in the Kerr metric. We show that the 
Schwarzschild/Kerr gravitational energy extracted over 
different regions of the space-time reproduces the ob¬ 
served bell-shaped behavior of the amplitude of the lower 
HF QPO [38M40L l42l] . The physical mechanism here pro¬ 
posed suggests that the behavior may be an effect related 
to an innermost marginally stable orbit predicted by GR 
in strong field regime, below which no stable orbital mo¬ 
tion can exist [11]. Because in theory such an orbit can 
have a generic eccentricity e [54| . in a more general way 
we refer to it as to innermost stable bound orbit (ISBO) 


instead of ISCO. 

The manuscript is organized as follow. In Section ITTI we 
review the main features of the Schwarzschild potential 
and give new hints on a possible way to reveal a signature 
of ISBO. Section M and Section m provide constrains 
on the maximum size allowed by tides of orbiting clumps 
of solid-state and plasma-state matter, respectively. In 
Section|V]we estimate both the energy loaded by tides per 
periastron passage on an orbiting clump of plasma and 
the amount of gravitational energy that would be emitted 
after tidal disruption of the clump. Section fVTl shows our 
results in the Kerr metric. In Section [VIl] we comment on 
other factors that might affect a signature of ISBO and 
show ('Section rVII Bl) an interpretation accounting for the 
observed coherence of the lower HF QPO as function of 
the luminosity of the source. Section IVIIII summarizes 
the conclusions. 


II. SCHWARZSCHILD POTENTIAL SHAPE 
AND ISBO 


We are interested in estimating differences of potential 
energy released by matter spiraling between the orbits in 
a curved space-time. For self-consistency we review some 
of the well known features of the gravitational potential 
in the Schwarzschild metric. Most importantly, we high¬ 
light new hints about a way of revealing a signature of 
ISBO. 

The familiar radial effective Schwarzschild potential 
per unit mass p. of an orbiting particle reads |54j 


v eff = 1 ~ 


2m 


2 mL 2 

[3 


P 

r 2 


(1) 


where m is the mass of the compact object and L is the 
angular momentunQ of the particle (per unit mass p) on 
the orbit of semi-latus rectum p and eccentricity e 


L = 


/ p 2 m 2 \ 1 ^ 2 

— 3 — e 2 ) 


( 2 ) 


The quantity p is defined such that the periastron of 
the orbit r „ = pm/(l + e) and for bound stable orbits 
p > 6 + 2e [541 ] . The second and third term in ([I]) refer 
to the gravitational potential in the Schwarzschild met¬ 
ric, while the forth term is the centrifugal potential. The 
attractive term oc 1/r 3 is relevant close to the compact 
object, i.e. when r ~ r g and causes the precession of 
elliptical orbits ib5i . 


1 The quantities here reported are expressed in geometric units 
G = c = 1, unless differently specified. The mass m in geometric 
units is equal to the gravitational radius of the compact object 
r g = GM/c 2 , where M is the mass of the compact object in SI 
units, G the gravitational constant and c the speed of light. For 
a 2 Mq neutron star r g ~ 3 km. 
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FIG. 1. Effective Schwarzschild potential as function of the 
radial coordinate r for two different orbits of given specific 
angular momenta L (blue curves). The periastra of the orbits 
are r p = 5.6 r g (ISBO) and r v = 8 r g . The apastra are at 
r a = 6.8 r g and r a = 9.8 r g , respectively. Red lines show the 
orbital energy E^ rb . 


FIG. 2. Difference of potential energy at periastron between 
two orbits with periastra 0.1 r g away. The energy expected to 
be released draws a bell-like shape, decreasing towards ISBO 
(r p = 5.6 r g ) and displaying its maximum in the range of radii 
r p ~ 6.5 — 10 r g . 


Fig. EO shows the known shape of the Schwarzschild 
potential for two different orbits. The blue lines draw 
the potential well for orbits of eccentricity e = 0.1 and 
periastra r p = 8 r g (top blue curve) and the ISBO r p = 
5.6 r g (bottom blue curve)@. The square of the orbital 
energy E 2 r u (in units of mass p) is reported (red lines) 
and reads [54j 


E 2 orb = 


(p — 2 — 2e) (p — 2 + 2e) 
p (p — 3 — e 2 ) 


(3) 


The abscissas of the crossing points between E° rb and 
V e ff are the periastron r p and the apastron r a of the orbit 
and come from the solution of the equation (dr/dr) 2 = 


Kb-Veff =0 


The ISBO at r p = 5.6 r g has a less 


deep minimum than the orbit at r p = 8 r g . It is because 
the term oc 1/r 3 in m starts dominating. 

The flattening of the minimum of the potential is a 
well known feature of GR [56]. Besides the precession of 
the orbits that has already been measured in the case of 
Mercury (but in a weak field regime), the observational 
implications of the term oc 1/r 3 in (jT]| in a strong field 
regime could be unmasked by the potential energy re¬ 
leased by spiraling matter close the compact object. It is 
well known from the accretion theory [57] that if matter 
accretes from an orbit at r = a r g and ends up onto an 
inner orbits at r = b r g (b < a), the energy released is 
the difference of potential energy between the two orbits. 
Such spiraling motion may be caused because of, e.g., the 
removal of orbital energy by tides acting on the clump of 
matter [58|. 

Fig. [2] shows a possible and simple way that could dis¬ 


2 For an orbit with e = 0 (ISCO) r = 6 r g . 


close a signature of ISBO. It depicts the behavior of the 
difference of the Schwarzschild potential, taken at perias¬ 
tron, between two orbits with periastra, e.g., 0.1 r g away. 
The orbits have a low eccentricity e = 0.1. The range of 
periastra is from r p = 30 r g to the ISBO at r p = 5.6 r g . 
Such difference of potential energy between the two or¬ 
bits, expected to be released, draws a bell-like shape: it 
increases up to ~ 7.5 r g and then decreases up to ISBO. 
This characteristic shape is a consequence of the term 
oc 1/r 3 in dTJ. In a Newtonian potential the differences 
of potential keep increasing in absolute value towards the 
inner regions. Thus, an observational evidence of the 
term oc 1/r 3 in |T]) and, therefore of ISBO, might be 
through these energetic considerations. We stress that 
these arguments have not yet been highlighted in the lit¬ 
erature. 

In Refs. [H 1U E3 it is shown that the energy (the 
amplitude) carried by the lower HF QPO in NS LMXBs, 
as function of its central frequency, draws a bell-shaped 
behavior like that around the maximum of fig. [2] The 
central frequency of the lower HF QPO is typical of the 
orbital motion within r ~ 6.5 — 10 r g . It may be worth¬ 
while to investigate whether the observed behavior orig¬ 
inates from differences of potential energy released by 
spiraling matter. 


III. SIZE ALLOWED BY THE STRONG TIDES 
OF ORBITING CLUMPS OF MATTER: THE 
SOLID-STATE CASE 

The arguments in Section [TT] give general clues about 
the overall pattern of the potential energy (in units of 
rest-mass energy) released by spiraling matter between 
orbits 0.1 r g away over different regions of the space- 
time, to compare it to the observed pattern reported in 
fig. 1 in Ref. [4^], fig. 3 in Ref. [38] and fig. 2 in Ref. [iof . 
We investigate whether this mechanism could actually 
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account for the observed energy carried by the lower HF 
QPO (10 34 —10 36 erg/s, from Atoll to Z-sources |4Cj). We 
need to give an estimation on the mass of the clumps of 
matter. The strong tidal force around a compact object 
sets an upper limit on their size. 

Consider a spherical clump of matter of radius R and 
its cup of height h = R /10 at the distance ~ R from 
the center of the clump. The volume of the cup is V = 
7 T h 2 (R — /i/3), thus its mass is p! = pV, where p is the 
density of the material. The effective gravitational force 
acting on the cup at distance r — R from the compact 
object is 


. 9 / 2m 
= pc* 1 


(r~R) 


6mL 2 


(r - R) 2 (r - R) 4 


(4) 


2 L 2 \ 

(r — R) 3 J 


where V e f / is from {T} and c 2 has been inserted to convert 
from geometric to SI units. The same cup placed on the 
opposite side of the clump at r + R would undergo the 
force F e ff = p'c 2 ( V ?* f ) , therefore the tidal force 

over the clump of matter is (see also Ref. [591] ) 


771 / 2 

ft = fi c 



(r-R) 


« p'c 2 2R 


d 2 Veff 
dr 2 


r 



(r+R) 


(5) 


where « indicates that we used a Taylor sum in R./r up 
to the first order since we expect R « r. The tidal 
force can not be larger than the ultimate tensile strength 
a of the material! (times the area 27r Rh of the cup), i.e. 
Ft < 2tt Rha. From this inequality we can get some 
order of magnitude on the maximum radius R allowed 
by tides 



( 6 ) 


where we write the density p = Y/c 2 , Y is the Young’s 
modulus of the material and c s the speed of sound in 
it. In the case of a solid-state material, e.g. asteroidfl 


3 In our solar system there are moons orbiting the parent planet 
below their gravitational Roche radius. They are held together 
by electrochemical bonds stronger than their self-gravity. Metis 
and Pan are moons orbiting below their gravitational Roche ra¬ 
dius. 

4 Even if an asteroid would probably be vaporized by the strong 
X-ray field (T ~ 10 7 K) far away the compact object, we consider 
this case for the sake of completeness. 


iron/steel-like material has a typical ratio in the labora- 
tory a/Y ~ 5 x 1CT 3 (cr ~ 10 2 MPa, Y ~ 10 2 GPa [H) 
and c s = 5 x 10 5 cm/s. From © the upper limit on 
the radius R of a solid-state clump that can survive to 
tides, at a distance r ~ 6 — 10 r„ from a 2 Mq neutron 
star, is ~ 30 cm (see also Ref. [6l[). The rest-mass en¬ 
ergy of a such sphere made of iron or steel (p ~ 8 g/cm 3 ) 
is pc 2 ~ 10 27 erg, thus from fig. [2] only ~ 6 x 10 23 erg 
can be emitted at the maximum of the bell-shape curve. 
This is too low for the observed energy of the lower HF 
QPO whose amplitude is some percent of the luminosity 
of the source, i.e. 10 34 — 10 36 erg/s. Hence, solid-state 
clumps made of iron/steel-like material (e.g. asteroids) 
are disfavored as sources of the lower HF QPO. 


IV. SIZE ALLOWED BY THE STRONG TIDES 

OF ORBITING CLUMPS OF MATTER: THE 
PLASMA CASE 

The plasma material in the accretion disk has different 
mechanical properties than solid-state matter. However, 
it may raise some doubts thinking that clumps of plasma 
in the accretion disk can survive to the strong shear forces 
due to differential rotation, thus not able to produce the 
lower HF QPO. On the other hand, to our knowledge, 
detailed simulations showing that clumps of plasma can 
not produce QPOs in the power spectra are not reported 
in the literature. Instead, there is a numerical code [52j| 
showing that orbiting clumps around a compact object 
can produce power spectra with twin-peak QPOs [531 ]. 
In this numerical code the clump is a collection of free 
particles. Each particle has its own geodesic. By means 
of ray-tracing techniques in the Schwarzschild metric the 
code shows that the signal from the clump, disrupted 
by both shear forces and tides, can produce twin QPOs. 
Most recently, Germana [(32| argued that the timing law 
of the azimuth phase <j)(t) on elliptical relativistic orbits 
could be the root for the appearance of multiple peaks in 
the power spectrum. Bakala et al. (jj modeled in detail 
the signal produced by orbiting hot-spots. Although the 
issue of a lower HF QPO produced by orbiting clumps in 
the disk may be of open debate, it is a viable interpreta¬ 
tion. We should bear in mind that the observed central 
frequency of the lower HF QPO is typical of the orbiting 
matter close to the compact object ITItI. [2~i1| . Moreover, 
the numerical code in Ref. [52| shows that already clumps 
of free particles can produce a lower HF QPO [53]. In 
reality, since the clumps of plasma in the disk might have 
some internal force, it is even more probable that they 
produce the lower HF QPO. Thus, it might be worth¬ 
while to investigate the mechanical properties of clumps 
of plasma in the accretion disk. 

We estimate the density of the inner part of the ac¬ 
cretion disk by using the equation of the density of the 
disk reported in Ref. [57]. For typical accretion rates in 
NS LMXBs of ~ 10 18 g/s (equation 1.7 in Ref. [HjJ; a Z 
source with luminosity L = LEdd ~ 2.5 x 10 38 erg/s [ioj . 
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where Lsdci is the Eddington limit), for a 2 Mq neutron 
star and at a distance of r ~ 6 — 10 r g , the density of 
the material in the disk is ~ 10 g/cm 3 , like that of iron. 
Because it is expected to be a highly ionized plasma, the 
speed of sound c s is not like that in iron materials. We 
write |57j 


Cs = 


/ 'yZkT \ 1/2 

V mi ) 


(7) 


where 7 ~ 5/3 is the adiabatic index, Z the charge state 
(Z = 1 for a hot plasma), rrii the ion hydrogen mass, k 
the Boltzmann’s constant. For typical temperatures in 
the inner part of the accretion disk ofT ~ 4 x 10' Kit 
turns c s ~ 8 x 10 7 cm/s. 

From © the upper limit by tides on the size of the clump 
depends also on the ratio cr/Y. This quantity describes 
the malleability of the body to external deformations, 
or, following the Hooke’s law (cr = Ye), it tells us how 
much deformable the material is before breaking. In the 
following we explore the idea of characterizing clumps 
of plasma through the ratio cr/Y. This way of proceed¬ 
ing might be justified by the results reported by Corral- 
Santana et al. [63^ who discovered large structures in the 
accretion disk of an X-ray binary system. Such structures 
propagate through the disk. Thus, it is hard to think 
to a completely smoothed accretion disk, but rather it 
might be characterized by inhomogeneities propagating 
throughout it. If inhomogeneities of plasma exist, as re¬ 
ported in Ref. [63), they have to be held together by some 
internal force (e.g. Ref. ]64j) against tides and differen¬ 
tial rotation. Whatever the nature of this force is, its 
value per unit area displays an ultimate tensile strength 
of the material. Equivalently, such force can be charac¬ 
terized through the mechanical binding energy, defined 
as the energy required to disaggregate/separate a body 
[651 . The binding energy per unit volume is the ultimate 
tensile strength a [ 66 ]. 

In this context, it might be worth noting that some 
features of the accretion disk theory is still of open de¬ 
bate. King et al. [13, [13 pointed out the yet unclear un¬ 
derstanding about the viscosity in accretion disk models. 
They investigated the disagreements between the viscos¬ 
ity parameter a required from observations and the one 
from simulations. King et al. 67] remark that a is not 
unambiguously the viscosity in the Navier- Stokes equa¬ 
tion but rather we can think it as the ratio of the stress to 
the rate of strain of the plasma. Its real physical meaning 
remains unclear. 


A. Constrain the ratio cr/Y of clumps of plasma in 
the accretion disk 

We estimate the ratio a/Y recalling the binding energy 
of a body. The gravitational binding energy of an uniform 
sphere of mass M and radius R is defined as the energy 
required to separate it in small parts and bring them at 


an infinite distance. The formula reads Et, = 3GM 2 /5R 
[65j |. This is valid for a body of loose material held to¬ 
gether by only gravity. We are dealing with clumps of 
hot-dense plasma, held together by some internal force. 
If cr is the ultimate tensile strength to be applied to break 
a specimen of such material, i.e. an energy per unit vol¬ 
ume [13, the energy required to disaggregate a body of 
volume V is the mechanical binding energy 


Eb = aV = -TrR 3 a 

O 


( 8 ) 


Substituting R from © we obtain 


cr 

Y 


(_3_ (J_) 3/2 f 1 _±) 3/2 ^b 

y 47 t \ 10/ ^ 30y Y 



(9) 
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Accreting clumps of plasma once disrupted by tides 
should release an energy that is at least of the same order 
of Eb- Thus, such binding energy may be of the order of 
that emitted by the lower HF QPO. For a Z source, the 
amplitude of the lower HF QPO is ~ 2 — 4% the luminos¬ 
ity of the source [4C| . We shall write Eb ^ 5 x 10 36 erg. 
Knowing the Young’s modulus Y = pc 2 s for the clump of 
plasma, we get cr/Y ~ 70. A ratio cr/Y > 1 means a 
deformation at break as big as more than once the initial 
size of the body. This is the case for, e.g., an elastic rope. 
Thus, a cr/Y ~ 70 means that the clump of plasma is far 
more deformable than solid-state iron materials. Because 
of this higher malleability to external deformations, © 
gives a bigger upper limit on the size of the structure of 
plasma that sustains tides, R ~ 6500 m. 

In Sec. m we mentioned that already clumps of free 
particles orbiting a compact object can produce the lower 
HF QPO [53j . Here the clump of plasma is held together 
by an internal force characterized by cr. We calculate 
whether the large upper limit of R ~ 6500 m set by tides 
can sustain the strong shear forces due to differential ro¬ 
tation. We divide the structure into, say, 10 slices and 
calculate the differential velocity as 


dv = 2 tt r 


( GM_\ 
\ 47 r 2 r 3 ) 


1/2 


2nr' 


/ GM \ 1/2 
\,4ttV 3 ) 


( 10 ) 


where r is the radius of the orbit at the center of the 
clump, r' = r -f- dh is the orbital radius of a contiguous 
slice. The shear force between two contiguous slices is 

F = pukdVdv (11) 

where p ~ 10 g/cm 3 is the density of the plasma, dV = 
2rrR 2 dh is the volume of the slice and 17 , the Keplerian 
frequency. It turns out that such load is F = 10 24 N. 
The ultimate tensile strength a ~ 70Y gives an opposite 
force keeping the slices together of 2ttR 2 <j ~ 10 25 N. 







6 


V. ENERGY LOADED BY TIDES ON 
ORBITING CLUMPS OF PLASMA 


With the limits on the size of the clump of plasma es¬ 
timated above, we may now give constrains on both the 
energy that tides might deposit on the orbiting clump 
and the amount of removed orbital energy m per pe- 
riastron passage. This will give us information on the 
interval of orbits the clump spirals over before disrupting 
and constrain how much potential energy is released, to 
compare it to the observations. 

The rate of energy deposited by tides on the orbiting 
clump of plasma is [69j 

= / dx 3 p v ■ VV eff (12) 


where v is the velocity of an element of the clump because 
of the perturbation by tides. We shall write |v| = dR/dt , 
with R radius of the clump. We turn the integral in time 
into an integral over the orbit, multiplying by dt/dr 


dE 

dr 


f „ dR 
= / dx P 


(13) 


The term dR/dr simulates the deformation of the clump 
because of tides as function of the distance r from the cen¬ 
tral object. On a slightly eccentric orbit dR/dr oscillates 
because of the changing tidal force between periastron 
and apastron. At r ~ 8 r g the tidal force © per unit 
area over a clump of plasma R = 3200 m is <7t ~ 10 17 
Pa, while the ultimate tensile strength of the material 
a = 70Y ~ 5 x 10 17 Pa. Thus, the clump may undergo 
some elastic phase before breaking- The tidal stress <tt 
slightly oscillates (by a factor ~ 2) from periastron to 
apastron on a time-scale @ t = \/v r ~3x 10~ 3 s. 
Because gt < a to get order of magnitudes we treat 
the clump as simple elastic body deformed by tides. We 
make use of the Hooke’s law da = Y dR/R , i.e. an in¬ 
finitesimal variation of the load causes an infinitesimal 
relative deformation proportional to the Young’s modu¬ 
lus Y of the material. The term dR/dr in (fl3l) turns into 
( R/Y)da/dr and since the load a is the tidal force © 
per unit area 


dE _ 1 

dr 4nR Y 


dx 3 p 


dFx 

dr 


VV e 


eff 


(14) 


The tidal force is constant over the clump at a given 
r. The integral over the volume, after transforming it in 



FIG. 3. Energy loaded by tides on a clump of plasma of radius 
R = 3200 m as function of the radial phase x of the orbit. 


spherical coordinates, reads « 4nR 3 (dV e ff /dr) r in the 
approximation R « r. Substituting © into m and 
integrating over the radius of the orbit, we obtain an 
estimation on the energy loaded by tides on the clump 


w , 3 tt pc 4 R 5 6 

~ 6.6 x lO^r 9 X 

(m 4 (21wr - 35m 2 - 3r 2 ) - 42mV+ 
— 15L 2 mr 2 (26m — 7r)^ 

We insert into (fl5l) the parametrized radius [54] 


(15) 


r (x) 


pm 

1 + e cos (x) 


(16) 


Fig.[3]shows the behavior of the energy loaded by tidet0 
(USD on the clump of plasma as function of the radial 
phase x- The curve is symmetric with respect to the 
apastron at x = 7T - The orbit has the periastron at 8 r g 
and eccentricity e = 0.1, around a 2 M© neutron stai@. 
Such energy is transferred from orbit to internal energy 
through tides [58, [69]. A removal of orbital energy of 
~ 10 36 erg from © gives a spiral motion of 0.4 r g ~ 1200 
m after the first periastron passage on these orbits with 
low e = 0.1. The tidal wave propagates through the 
clump on a time-scale r = 2 R/c s . The radius of the 
fraction of the clump disrupted between two periastron 
passages is Ri = tc s /2 ~ 1100 m. On the second peri¬ 
astron passage the radius of the clump might be smaller 
(R ~ 2000 m), the energy deposited by tides (fl5l) now is 
E(r) ~ 10 35 erg causing a spiral motion of 0.1 r g ~ 300 
m. On the third periastron passage what is left has a 


5 Structures equal/larger than the limits constrained in Sec IIV Al 
can not survive at all because of tides. On the other hand, struc¬ 
tures smaller than R ~ 3000 m would not be able to emit the 
required energy seen in the lower HF QPO in Z sources, i.e. 

5 x 10 36 erg/s. 

6 is r is the relativistic radial frequency of the orbit. In the 
Schwarzschild metric around a 2 M® neutron star u r ~ 350 Hz 
at its maximum (r ~ 8 r 5 ). 


7 Note that the order of magnitude estimated here 
E(r) ~ 0.12% fie 2 agrees with that from the formalism de¬ 
scribed in Ref. [4^1 for a star disrupted by a supermassive black 
hole. 

8 All over the window of periastra ~ 5.6 — 10 r g the energy de¬ 
posited on the clump R = 3200 m is ~ 10 36 erg. 
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FIG. 4. Difference of potential energy at periastron between 
two orbits with periastra ~ 1700 m away as function of the 
orbital frequency. To compare it to the data see fig. 1 in 
Ref. [42], fig. 3 in Ref. [38| and fig. 2 in Ref. [40 l], 


radius R ~ 1000 m < Ri. The energy deposited by tides 
on such relic at each periastron passage is E{r) ~ 10 33 
erg causing a spiral motion of 0.01 r g ~ 30 m. The bind¬ 
ing energy © is Eb ~ 10 34 erg. Thus, the number of 
periastron passages before disrupting may be ~ 10. In 
total the clump might make ~ 12 periastron passages. 
The number of Keplerian turns if] ~ 24. The difference 
of potential energy between the orbits is emitted by such 
an event on a time-scale r' =12/(350 Hz)~ 0.034 s. The 
coherence as defined in Section |T| Q = v/Av in this case 
is Q = v k r' ~ 24, corresponding to the number of Kep¬ 
lerian turns. 

In total, the initial clump with R = 3200 m has spiraled 
over an interval of orbits ~ 0.6 r g ~ 1700 m. Fig. [4] 
shows the differences of potential energy, at periastron, 
released between two orbits with periastra ~ 1700 m 
away as function of the orbital freauencTl around a 
2 Mq Schwarzschild compact object (the curve spans 
over 7' ~ 10 — 6 r g ). It displays the characteristic bell- 
shape centered at ~ 700 — 800 seen in the obser¬ 
vations [3814+1 [42| . The emitted energy matches that 
carried by the lower HF QPO in Z sources [40] and cor¬ 
responds to ~ 0.5% lie 2 . The figure corresponds to the 
maximum of the overall bell-shape in fig. [2] Because the 
bell-shape in fig. [4] is typical of the observed amplitude 
behavior of the lower HF QPO rather than the upper 
one, this study may add one more clue: it suggests that 


9 Around the region r~8r s the ratio of the relativistic Keplerian 
frequency to the radial one u r is v^/vr ~ 2. Thus, if the 
number of periastron passages is n, the number of Keplerian 
turns is 2 n. 

10 We calculated the orbital frequency for an orbit with e = 
0.1 by means of the integrals reported in Ref. [5^1 describing 
both the azimuth phase 0(x) and coordinate time t(x) in the 
Schwarzschild metric, for generic orbits of eccentricit y e . We 
resolved them by applying a Taylor sum around e = 0 [6^1 . 

11 For a 1.8 Mq compact object the curve in fig. [4] is centered 
around 800 Hz. 


the lower peak of the observed twin-peaks corresponds 
to the Keplerian frequency. This is in agreement with 
the results from the numerical code on tidal disruption 
(53|, where the lower HF QPO corresponds to the Ke¬ 
plerian modulation while the upper HF QPO to the 
modulation at + v r . Such conclusion is also reported 
in Ref. [H]. Models proposing other mechanisms to pro¬ 
duce HF QPOs and assuming the lower HF QPO to be v\~ 
have been previously proposed [l8l-[20j|. differently than 
in Refs. S S EE Ea in which the upper HF QPO is 
linked to v k . 

The lower HF QPO is observed to range ~ 600 — 950 
Hz mm IH . From fig. [4] we deduce that the clustering 
might be because the difference of potential energy re¬ 
leased by such random events between near orbits has its 
maximum in the region of the space-time r ~ 10 — 6 r g 
(see also fig. The steeper decrease in fig. [4] is drawn 
by the approach to ISBO (see also fig. m- Hence, the ob¬ 
served behavior of the amplitude of the lower HF QPO as 
function of its central frequency may give clues related 
to the ISBO predicted by GR in a strong field regime 
(r ~ r g ). 


VI. THE KERR METRIC CASE 


We performed our calculations in the case of a rotating 
compact object, whose space-time around it is described 
by the Kerr metric [HI. The effective radial gravitational 
potential V e ff t Kerr is such 0 that for a non-rotating 
compact object (a = 0, a is the angular momentum per 
unit mass of the compact object) we get JTJ 


V eff, 


Kerr 


= l - 


2m 


a 


4 - 77 + 


+ 


r r‘ 
4amE or bL 


1 - 


2 m\ a 2 E 2 rb 
r ) r 2 
2 m\ L 2 


+(17) 


/ 2 

r J r z 


with 0 < a < m and E or b , L orbital energy and specific 
angular momentum of the orbiti ng p article on a generic 
circular orbit in the Kerr metric [131 ] . 

We note that the Kerr parameter a lowers the maxi¬ 
mum radius R of the clump of plasma allowed by tides: 
for an orbit at 8 r g R ~ 7000 m [3 for a = 0, R ~ 6500 
m for a = 0.9m. The work done by tides per perias¬ 
tron passage remains of the same order of that obtained 
in the Schwarzschild metric (equation (fThl) . E(r) ~ 10 36 
erg). Thus, the clump with size R = 3200 m as in Sec. IVl 
makes an almost equal number of turns before disrupt¬ 
ing, spiraling over an interval of orbits of ~ 0.6 r g ~ 1700 
m. Fig. [5] shows the difference of Kerr potential energy 
released for different values of a between orbits ~ 1700 
m away. We see that such mechanism is more efficient 
for increasing a: the amplitude of the bell-shape curve 


12 This value agrees to that from © for an orbit with e = 0. 








v (Hz) 

FIG. 5. Differences of Kerr/Schwarzschild potential energy 
expected to be released by an orbiting clump of plasma of 
radius R = 3200 m after being disrupted by tides. The curves 
refer to different values of the Kerr parameter a. 

increases. The effect of a is also to move the ISBO to 
the inner regions, so that the orbital frequency extends 
to higher values. The figure refers to a 2 M 0 compact 
object. For a 2.5 M 0 each curve in the figure is almost 
shifted to lower frequencies by ~ 200 Hz. It is known 
that different couples (to, a) could describe the same fre¬ 
quency at a given orbital radius [32j. Instead the value of 
the mass has no effects on the amplitude of the bell-shape 
curves in fig. [3] 


VII. OTHER FACTORS AFFECTING THE 
AMPLITUDE OF THE LOWER HF QPO 

We advise caution about the behavior_shown in figs. [4] 
[5] and seen in the observations [38l - l40l . [42j as candidate 
to disclose a signature of ISBO. We discuss below other 
factors that may affect the amplitude of the lower HF 
QPO. 


A. The dependence of the amplitude on the energy 

The amplitude of the lower HF QPO measured at dif¬ 
ferent energy channels keeps increasing towards hard X- 
ray. At ~ 20 — 30 keV the amplitude of the lower HF 
QPOs is ~ 20%, while at ~ 2 keV is less than 5% [42l.[70j. 
The lower HF QPO is usually stronger in the hard state 
of the source, where the energy spectrum is dominated by 
a power law. This means that the thermal emission from 
a Keplerian disk can not justify its nature. A corona of 
hot electrons scattering to higher ener gy t he seed pho¬ 
tons from the disk has been proposed [7ll - l74l |. On the 
other hand, Schnittman [75[ noted that this mechanism 
may not amplify the modulation, but rather would lower 
the amplitude of the modulation at high energies, since 
the scattering of photons by hot-electron would smooth 
the oscillation from the disk. 

Alternatively to the corona scenario, the modeling in 


Ref. [7(| shows that sub-Keplerian motion in the inner 
part of the accretion disk (r < 20 r g ) might cause a 
postshock region with hard X-ray emission. Hard X-ray 
emission also might come from the converging bulk mo¬ 
tion of the accreting matter. In Ref. [77} it is shown that 
in such sub-Keplerian region the specific angular momen¬ 
tum of the accretion flow can also behave in a Keplerian 
way. Therefore, Keplerian imprints like that described in 
our manuscript may be seen in the observed hard state. 
Yan and Wang [73] pointed out the idea that the hard 
X-ray spectrum of the lower HF QPO might originate 
by synchrotron mechanisms from magnetized hot-spot. 
Calculations on the radiation that would be emitted by 
a tidally disrupted magnetized clump of matter were at¬ 
tempted in Refs. [521 l79f ] and full magnetohydrodynamic 
simulations are demanded. Thus, the exact energy emis¬ 
sion mechanism of the lower HF QPO may remain of 
open debate and is far beyond the goal of this paper. 
This manuscript attempts to propose a mechanism jus¬ 
tifying how the quantity of energy carried by the lower 
HF QPO would originate and not the energy spectrum of 
the emission, which depends on the exact nature of the 
energy emission mechanism. We are comparing our cal¬ 
culations to bolometric luminosity. Gravitational energy 
is extracted through tides and transformed into radiation 
through a given emission mechanism. The bolometric lu¬ 
minosity emitted should be equal to the amount of gravi¬ 
tational energy extracted. Thus, the behavior in figs. EE] 
and seen in the observations [381 - 4401 1421 ] may still remain 
a candidate to disclose a signature of ISBO. 


B. The dependence of the amplitude on the 
luminosity of the source 

It is known that the amplitude of the lower HF QPO 
is different among sources with different luminosity, i.e. 
accretion rate. The amplitude of the lower HF QPO is 
up to ~ 20% in atoll sources, while in Z sources is ~ 5%. 
Z sources have a luminosity (~ 10 38 erg/s) that can be 
up to two orders of magnitude higher than atoll (~ 10 36 
erg/s) [10]. The unique source XTE J1701-462 displays 
both phases, i.e. atoll and Z [80]. In fig. 4 of Ref. [[80] is 
shown the amplitude (in percent of the total luminosity 
of the source) in the atoll and Z phase. The same source 
displays different amplitudes of the lower HF QPO (10% 
in the atoll phase, 3% in the Z phase) at the same central 
frequency of the peak. Thus, it has been suggested that 
the amplitude of the lower HF QPO as function of its 
central frequency within a source (e.g. fi gs. E El in this 
paper, fig. 1 in Ref. [42], fig. 3 in Ref. |38| and fig. 2 
in Ref. [40j) may not be simply used as clue related to 
a signature of ISBO [80|. The study proposed in this 
manuscript may offer a possible interpretation to these 
differences between atoll and Z sources. 

In the Z phase XTE J1701-462 has a luminosity of 
L ~ 0.5 LEdd. ~ 10 38 erg/s (fig. 5 in Ref. [80|l. where 
LEdd ~ 2.5 x 10 38 erg/s is the Eddington luminosity for 
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FIG. 6. Simulated quality factor Q (coherence) of the Keple- 
rian modulation by an orbiting clump of plasma as function 
of the luminosity L of the source (in units of the Eddington 
luminosity Lsdd ~ 2.5 x 10 38 erg/s) for two representative 
sources (see text). Such coherences are typical of the lower 
HF QPO in these sources (see fig. 5 in Ref. [80| and fig. 3 in 
Ref. |40[ for a comparison with the data). 

a ~ 2 Mq neutron star [57]]. The accretion rate of the 
source is M ~ 7 x 10 17 g/s. The lower HF QPO has 
a coherence up to Q ~ 30 and amplitude ~ 3%, corre¬ 
sponding to ~ 3 x 10 36 erg/s. At lower accretion rates 
(M ~ 7 x 10 16 g/s) in the atoll phase the coherence is 
up to Q ~ 150 and the amplitude up to ~ 10%, also 
corresponding _to ~ 2 x 10 36 erg/s (for such details see 
fig. 5 in Ref. 0). From the calculations in Sec. IIV Al 
and Sec. E we obtain that in the Z phase structures as 
large as R ~ 3000 m can survive to tides, making 10 peri- 
astron passages, producing an oscillation with coherence 
Q ~ 20 (typical of the Q observed in XTE J1701-462 
in its Z phase 0 ) and emitting the observed energy 
of the lower HF QPO ~ 3 x 10 36 erg, corresponding to 
~ 0.6% pc 2 . The limiting size Ri given by the accretion 
rate is Ri ~ 900 m. The initial clump with R = 3000 m 
spirals over an interval of orbits of 0.73 r g ~ 2100 nr. 

In the atoll phase of XTE J1701-462 the accretion rate 
is M ~ 7 x 10 16 g/s and the lower HF QPO has also 
a luminosity ~ 10 36 erg/s (10% of the luminosity of the 
source). Because of the lower accretion rate, the density 
is lower by a factor ~ 5 (p ~ 1 g/cm 3 ) than in Z phase. 
A structure with R ~ 3000 nr can spiral over 1.5 r s and 
emits ~ 2 x 10 36 erg, corresponding to 1.3% pc 2 . The 
limiting size is Ri ~ 500 nr and this last stage lasts for 
~ 60 periastron passages, spiraling over 0.2 r g . Thus, the 
lower HF QPO has a coherence of Q ~ 130, much higher 
than that in the Z phase, as in the observations [80| . 

We did the calculations also for the source 4U 0614+09, 
which has the lowest luminosity in fig. 5 in Ref. [0 (see 
also fig. 3 in Ref. 0)- Its accretion rate isM~7xl0 15 
g/s and the maximum luminosity of the lower HF QPO 
is ~ 10 35 erg/s (~ 10% that of the source). The limiting 
size given by the typical speed of sound c s is Ri ~ 400 


m, the density of the plasma p ~ 0.5 g/cm 3 (one order 
of magnitude lower than Z sources). An initial structure 
with R ~ 2500 m spirals by 0.8 r g , 0.4 r g , 0.2 r g , 0.1 r g 
and 0.03 r g after the first, second, etc periastron passage. 
Afterwards the relic reaches the limiting size R ~ Ri 
and spirals by 0.003 r g per periastron passage making 
~ 15 passages. In total, it might make ~ 20 periastron 
passages before disrupting (Q ~ 40). The energy emitted 
is ~ 3 x 10 35 erg corresponding to ~ 1.2% pc 2 . 

In Z sources the clump makes a similar number of turns 
before disrupting. The allowed number of turns by tides 
might be dictated by the bigger size Ri , allowing tides to 
make more work per periastron passage. 

Fig. [6] shows the simulated coherence Q of the Keple- 
rian oscillation obtained with such mechanism for both 
the atoll and Z phase of XTE J1701-462 and the lowest 
luminosity source 4U 0614+09. The figure reproduces 
the observed coherence of the lower HF QPO as function 
of the luminosity (accretion rate) of the source (fig. 5 in 
Ref. 0 and fig. 3 in Ref. 0 ). These results suggest 
that the lower HF QPO might correspond to the Keple- 
rian modulation of orbiting matter around the compact 
object, similar to what the results in Sec. E suggest. 

From the mechanism proposed here, we see that struc¬ 
tures of plasma with size i? ~ 3 x 10 3 m (at an orbital 
radius r ~ 8 r g ~ 24 km, R/r ~ 0.1) are required to 
produce the observed luminosity of the lower HF QPO 
among a sample of sources with different accretion rates. 
They would emit an energy that is ~ 10 36 erg in both 
atoll and Z sources (~ 10% and ~ 3% the total lumi¬ 
nosity, respectively) and in the lowest luminosity state 
(e.g. 4U 0614+09 with density of the plasma one order 
of magnitude lower) ~ 10 35 erg. They might produce 
Keplerian modulations with different coherences similar 
to the observations (fig. [5]), because of the interplay be¬ 
tween tides and the mechanical properties of the plasma. 
Speaking about orders of magnitude, we notice that the 
size of the structure required to emit the observed en¬ 
ergy is similar among a sample of sources spanning ~ 2 
orders of magnitude in luminosity. The geometry of the 
space-time because of the neutron star might play a rele¬ 
vant role among the sample of sources. Larger structures 
may not survive/form at all because of tides. This may 
justify why we do not observe an amplitude of the lower 
HF QPO of, e.g., 10 — 20% in Z sources (~ 10 37 erg/s), 
or higher than those observed in atoll sources. 

From m we derive a/Y ~ 70 in the Z state of 
XTE J1701-462 and also for the Z source in Sec. IIV Al 
a/Y ~ 300 in the atoll states, i.e. both the atoll phase 
of XTE J1701-462 and 4U 0614+09, even if they differ 
by one order of magnitude in luminosity. Thus, the atoll 
and Z states of LMXBs may characterize some properties 
of the plasma in the accretion disk (e.g. electrochemical 
bonds and/or magnetohydrodynamical). 

To conclude, different accretion flows between Z and 
atoll sources might change the mechanical properties of 
the clumps of plasma, such as p, c s , Y and a/Y, in a way 
that both the amplitudes and coherences of the lower HF 
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QPO also depends on the accretion flow [80']. However, 
within a source that does not change dramatically its 
accretion rate the behavior of the amplitude of the lower 
HF QPO with its central frequency (i.e. fig s. SI [5] in this 
paper, fig. 1 in Ref. [lU, fig. 3 in Ref. [38j ] and fig. 2 in 
Ref. [40]), after this analysis, may still remain a candidate 
for a signature of ISBO. 

VIII. REMARKS AND CONCLUSIONS 

Reis et al. [49j reported a QPO in the X-ray flux 
from the tidal disruption of a star by a supermassive 
black hole. Here we have shown that the amplitude of 
the lower HF QPO in NS LMXBs might display the 
Schwarzschild/Kerr potential energy released by spiral¬ 
ing clumps of plasma in the accretion disk. By means 
of first-approximation calculations this mechanism gives 
orders of magnitude typical of the observed energy car¬ 
ried by the lower HF QPO. We would stress that magne¬ 
tohydrodynamic modeling is certainly required to reach 
firmer conclusions. The estimations reported here are 
promising and are shown for the first time. 

The mechanical properties of the plasma (c s , Y, a/Y) 
in the accretion disk depend on its density and temper¬ 
ature and, therefore, on the accretion rate of the source. 
This study suggests that both the amplitude and coher¬ 
ence of the lower HF QPO also might depend on the in¬ 
terplay between the mechanical properties of the clumps 
(c s , Y, <j/Y) and the work done by tides on them. Such 
differences among sources with different accretion rate 
are reported in Ref. [ifl [gcj and were also discussed in 
Ref. gj. 

The main and new results of this study are: (i) the dif¬ 
ferences of Schwarzschild/Kerr potential energy emitted 
by spiraling clumps of plasma might account for the ob¬ 
served energy (amplitude) carried by the lower HF QPO. 
Also, the coherence time obtained from the time-scale of 


tidal disruption (~ 0.03 — 0.2s) is typical of the observed 
ones (~ 0.01 — 0.2 s; Refs. [38H401 [80] 1. This is the first 
time that a proposed physical mechanism can justify both 
the observed emitted energy and coherence of the lower 
HF QPO. We are able to give an interpretation to the dif¬ 
ferent coherence of the lower HF QPO among a sample 
of sources with different luminosity (fig.®- (ii) Within a 
source with a given accretion rate, the typical bell-shape 
of the amplitude of the lower HF QPO as function of its 
central frequenc y (s een in the data collected from sev¬ 
eral NS LMXBs [38^4(1 EH) might be reproduced by dif¬ 
ferences of Schwarzschild/Kerr potential energy between 
near orbits over the region r ~ 6 — 10 r g . Thus, such bell¬ 
shaped behavior may be related to the terms oc 1/r 1 2 3 4 5 6 7 in 
© and m in a strong field regime and therefore could 
be a geometry-related effect due to ISBO (figs. [2j ® and 
®. We however would stress that the Schwarzschild and 
Kerr metrics only approximate the space-time around a 
neutron star & (iii) Our calculations might justify for 
the first time the interval of radius where the lower HF 
QPO is produced, i.e. ~ 6—10 r g : In this region the high¬ 
est difference of potential energy is released between the 
orbits the clumps spiral over before disrupting (fig. [2J. 
Future observations with a higher signal to noise ratio 
[8^ ] could extend this region which, with the data to 
date, it may just be masked by the background radiation 
from the disk. 
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